NQR and X-ray investigation of the structure of Na 2 /3Co02 compound 
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We have synthesized various samples of the x = 2/3 phase of sodium cobaltate Na^CoC^ and 
performed X-ray powder diffractions spectra to compare the diffraction with the structure proposed 
previously from NMR/NQR experiments [H. Alloul est al, EPL 85, 47006 (2009)]. Rietveld analysis 
of the data are found in perfect agreement with those, and confirm the concentration x=2/3 obtained 
in the synthesis procedure. They even give indications on the atomic displacements of Na inside the 
unit cell. The detailed NQR data allow us to identify the NQR transitions and electric field gradient 
(EFG) parameters for 4 cobalt sites and 3 Na sites. The spin-lattice and spin-spin relaxation rates 
are found much smaller for the non-magnetic Co 3+ sites than for the magnetic sites on which the 
holes are delocalized. The atomic ordering of the Na layers is therefore at the source of this ordered 
distribution of cobalt charges. The method used here to resolve the Na ordering and the subsequent 
Co charge order can be used valuably for other concentrations of Na. 

PACS numbers: 71.27.+a, 61.66.f, 76.60.Gv 
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I. INTRODUCTION 

Since the discovery of high thermo-electric power 
(TEP) 2 and superconductivity^ in Na cobaltates, ex- 
tended efforts have been done in order to understand 
their magnetic and electronic properties. Anomalous 
magnetic properties have been discovered in unexpected 
ranges of Na concentrations, with an abrupt change of 
magnetic correlations from AF for x < 0.62 to ferromag- 
netic for x > 0.67—. However, a situation quite unusual 
in solid state physics prevailed so far, as very few experi- 
ments have allowed to correlate these magnetic properties 
with Na ordering, apart for the featureless band insulator 
for x = 1— £,and the very peculiar x = 1/2 phased - —. 

We have recently reported that a combination of nu- 
clear magnetic resonance (NMR) and nuclear quadrupo- 
lar resonance (NQR) experiments did concur to let us 
determine the Na order in the x = 2/3 phaseA, which 
is a nearly ferromagnetic metallic phase which does not 
order magnetically down to the lowest temperatures^. 
For this particular phase the large set of NMR data ob- 
tained in the past allowed us, when complemented by 
NQR data, to demonstrate finally that the Na order is 
rather simple. It results in a differentiation of Co sites 
into non-magnetic Co 3+ sites in which the t2 fl subshell of 
electronic levels is filled, and a metallic Kagome network 
of Co sites on which the doped holes are delocalizedA. 

We perform here a careful Rietveld analysis of the 
X-ray data which confirms the structure proposed from 
NMR and NQR and definitively establishes the x = 2/3 
Na content of this phase contrary to the x=0.71 value es- 
timated by others from chemical analyzes of single crys- 
tal samples^!. Furthermore, we report in detail here the 
NQR data which, together with the approach developed 
m Ref. [H will help in similar analyzes for other struc- 
turally more complicated cases detected for different x 
values. 

The paper is organized as follows. In section II we 
display the methods used to synthesize single phase sam- 



ples for x = 2/3 and display the powder X-ray spectra in 
which the Bragg peaks associated with the Na ordered 
substructure are resolved. The best Rietveld fits allowed 
us to determine the slight atomic displacements of the Na 
atoms in the model unit cell deduced from NMR/NQR 
data. In section III we report the Na and Co NQR spec- 
tra and show that we could detect most of the NQR lines 
associated with 4 Co sites identified by the former com- 
bination of NQR and NMR data. Some of the low fre- 
quency lines were however difficult to detect owing to ex- 
perimental limitations. At the end of section III we show 
that the accurate electric field gradient (EFG) parame- 
ters can be hardly described merely by a point charge 
model, even taking into account the site displacements 
determined by X-rays. This reveals the need for more 
sophisticated and accurate calculations taking fully into 
account the atomic and electronic band structure. The 
spin-spin and spin-lattice relaxation data reported in sec- 
tion IV for the four Co sites are studied and analyzed 
then, which allowed us to establish that both quantities 
can be used to distinguish the non-magnetic cobalts from 
those of the more magnetic Co sites which are organized 
in a Kagome sublattice in the x = 2/3 phase. 



II. SAMPLE PREPARATION AND LATTICE 
STRUCTURE 

It has been established that the phase diagram of Na 
cobaltates is discontinuous and consists of a series of ho- 
mogeneous phases extending over limited concentration 
ranges, separated by composition gaps^IIi. A homoge- 
neous domain range is characterized by a specific Na or- 
dering which one can monitor by powder X-ray diffrac- 
tion experiments. These diffraction patterns were taken 
at 300 K in the range of 26 angle 10-130° (Cu Ka radia- 
tions were used) . X-ray spectra such as that displayed in 
Fig. [1] allow to determine the lattice parameters from the 
X-ray diffraction peaks of the hexagonal structure of the 
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FIG. 1: (color online) Na 2 /3Co02 structure determination. Main field: X-ray diffraction pattern observed I b s in blue, calculated 
Icaic in red and I bs — Icaic in green. The intense 002 substructure reflection at 29 ~ 16° is not used for the refinement calculation 
because of its strong instrumental asymmetry. Insert: Intensities of satellite reflections: blue observed, red calculated. 



C0O2 substructure. Also, as shown in the inset of Fig. [TJ 
a set of satellite reflexions associated with the specific Na 
ordered structure could be detected and were found quite 
distinct for the various stable phases. 



A. Sample synthesis 

We have evidenced that in pure oxygen atmosphere the 
compounds with x — 2/3 are stable in the temperature 
range 600-700°C. Phases with higher Na content are sta- 
bilized at higher temperatures, which for a nominal com- 
position x — 2/3 releases some C03O4. Conversely lower 
preparation temperatures favor phases with x < 2/3, and 
in such case the excess Na is somewhat difficult to detect 
by diffraction, although after treatment in air the charac- 
teristic diffractions of Na2C03 appear. So to retain the 
x = 2/3 phase in powders, the samples were directly in- 
troduced in a furnace stabilized at 700° C, and were then 
quenched at room temperature after treatment. Three 
different approaches have been used to synthesize single 
phase materials: 

1) direct synthesis from a stoichiometric composition 
of C03O4 and Na2C03. However we found that the reac- 
tion is rather slow at such low temperatures and requires 



weeks of treatment with intermediate grindings to elimi- 
nate the last traces of unreacted C03O4. 

2) from a mixture of cobaltates with calibrated com- 
positions synthesized previously (such as Nai/2Co02 and 
Nao.7iCoOz^ taken in a proper ratio. 

3) by de-intercalation of Na from Nao.7iCo02 by an- 
nealing it at 700° C - out of its own stability temperature 
range. Of course, in that case the composition x = 2/3 
of the material prepared differs from its nominal compo- 
sition x — 0.71, which remains unchanged and would be 
detected by chemical analysis. The released Na remains 
in excess and is not well crystallized (oxidized) and diffi- 
cult to observe by X-rays. 

Overall the second method was found the most reliable 
and was easily reproducible. 



B. Rietveld refinement of the structure 

As already indicated in a previous published report^, 
the NMR and NQR data have allowed us to identify 3 
Na and 4 Co sites, their occupancy and those with axial 
symmetry in the lattice. It has been then possible to es- 
tablish the lattice cell of the structure displayed in Fig. [5] 
Such a structure ought to be recovered by X-ray scatter- 
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FIG. 2: Na 2 /3Co02 structure model. Two neighboring 
sodium and cobalt layers are shown in one unit cell size. The 
upper plane (a) illustrates the Na ordering and the lower one 

(b) the differentiation in the Co plane. The dotted circles rep- 
resent the displaced positions of the sodium ions, as deduced 
from the X-ray data analysis. The diameters of the Na and 
Co ions scale with the atomic radii of the corresponding ions. 

(c) Top view of the two considered Co and Na planes. 



ing experiments and should explain the set of satellite 
diffractions detected, such as those displayed in the inset 
of Fig. [TJ 

For the phases with x > 0.5 we have detected such 
sets of satellites, which in most cases can be described 
by a single incommensurate wave vector with component 
q%. Under these conditions the lattice loses the hexago- 
nal symmetry and becomes generally orthorombic with 
a or t = a-hex^S, b ort = a hex /q£, c ort = c hex , as was no- 
ticed early by Zandbergen et al. from high resolution 
microscopy experiments^. This has led to propose a phe- 
nomenological models in which an increase of Na content 
corresponded to the insertion of " stripes" of Na2 with in- 
creasing a ort in the known structure of Na^CoCv It is 
easy to evaluate that in such a model the empirical rela- 
tion x = 1 — ql applies. 

In contrast only the particular composition x = 2/3 ex- 
hibits a commensurate locking leading to a conventional 
superstructure with a rhombohedral unit cell a r h — 
^ a hex^/i, Crh = 3 c hex which, in the phenomenological 
description, corresponds to ql = 1/3, as expected. This 
structure, which corresponds to the spatial group R-3c 
(n°167) allowed excellent quality Rietveld fits of the X- 
ray diffraction pattern as can be seen in the expanded 
diagrams of Fig. QJ£. These Rietveld refinements lead 
to the distribution and occupancy of cobalt and sodium 
sites given in Table [H in perfect agreement with the 
NMR-NQR results. Moreover, they definitively confirm 
the composition x — 2/3. Indeed, only the sodium sites 
mentioned in the Table U are fully occupied whereas the 
other sites deduced from the hexagonal substructure re- 
veal themselves to be completely empty within the limits 
of standard deviations. 



TABLE I: X ray determination of the structure of Na 2 /3Co02 , 
data collection 26 = 10-130°, Cu K a . a =9. 8007(1) A; 
3c =32.8151(4)A; spatial group R-3c (n° 167), hexagonal 
axes, Z = 72, R = 5.41%, R w = 7.16%. 
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0.187(2) 
1/6 6 
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l/2 b 
l/6 b 
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Ox-1 
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0" 


0.0298(5) 


1.3(2) 


Ox-2 


36 f 


o b 


0.0290(4) 




Ox-3 


36 f 


l/6 b 
l/2 b 


0.0320(5) 




Ox-4 


36 f 


0.0295(5) 





"Number of sites and WyckofF notation for position. All sites are 
found to be fully occupied (see text). 

''Those coordinates have been left free but the positions deduced 
from the fits did not deviate significantly within experimental ac- 
curacy with respect to the mean unit cell positions. 



Most of the ionic coordinates are unchanged with re- 
spect to the substructure positions (see Table A} . The 
most displaced are the Na2 ions. The displacement of 
Na2 sites already were mentioned in the Ref. [l2lfl3l but 
it was attributed to the repulsion of randomly located 
neighboring Na ions, locally violating the hexagonal sym- 
metry. But as one can see on Fig. [5] the relaxation of 
position of these sites corresponds to a dilatation of the 
Na2 triangles, which agrees with the interionic repulsion 
in this highly charged region. As might be expected, the 
displacements of Na along the honeycomb lattice (ionic 
conduction channel) of the outer Na2b triangle are nearly 
twice larger than those of the inner Na2a triangle. 

In addition we would like to mention that we have 
calculated the diffraction spectra for eight distinct sam- 
ples and found in all cases that the Cola position is 
slightly shifted along the c-axis direction towards the Na2 
trimers, with the same offset » 0.06(1) A. The oxygen 
ions are only found slightly displaced with respect to the 
ideal unit cell, which points the compactness of the C0O2 
slabs between which are inserted the Na ions. 



III. NQR SPECTRA AND EFG PARAMETERS 

A nucleus with nuclear spin I > 1/2 has an electric 
quadrupole moment in addition to its nuclear magnetic 
moment. The nucleus interacts with the electronic envi- 
ronment not only through magnetic hypcrfine couplings 
due to its magnetic moment, but also through the interac- 
tion of its quadrupole moment with the local crystal elec- 
tric field gradient (EFG), either static or dynamic. The 
EFG arises from a non-symmetric distribution of electric 
charge around the nucleus. This charge can originate 
from non-bonding electrons, electrons in the bonds and 
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charges of neighboring atoms or ions. Therefore NQR is 
a sensitive tool for studying solids as it provides detailed 
information on the static and the dynamic properties of 
the structure on the scale of a few interatomic spacings. 
Thus, NQR may be regarded as a powerful tool for in- 
vestigating the local order in solids whereas the inter- 
pretation of the data on disordered (or complex-ordered) 
materials from usual scattering experiments, such as X- 
ray or neutron scattering, is complicated by the absence 
of the long-range order translational symmetry or large 
and complicated unit cell. 

Therefore to demonstrate these advantages of NQR to 
the non-NMR specialist we start this section by a brief 
NQR background introduction (A) and describe in (B) 
the basic features of the experimental techniques used. 
Next we show the 23 Na NQR spectrum (C) and compare 
it to the 23 Na NMR datail We present then (D) the 
59 Co NQR spectrum and show how fast and slow relaxing 
cobalt sites can be separated in the spectrum. Finally 
in (E) we show that the accurate EFG parameters can 
hardly be described merely by a point charge model, even 
taking into account the site displacements determined by 
X-rays. 



The change of tuning frequency of the NQR probe has 
been done with a variable capacitor. However the Q- 
factor of the probe (and the sensitivity of the spectrom- 
eter) varies considerably with frequency. This limits se- 
riously the accuracy of absolute intensity measurements 
of NQR signals. Relative intensity measurements could 
only be done accurately for lines occurring in a narrow 
frequency range. 

Due to the spin-spin relaxation process (which will be 
discussed in detail later) the intensity of the observed 
spin echo NQR signal decreases with increasing delay r 
between the RF pulses. It is usually impossible to re- 
duce this time to zero, as the receiver is overloaded by 
the RF pulse and only recovers after a minimum time 
tr. Furthermore, for short values of r the detected signal 
is heavily distorted by residual oscillations (called "ring- 
ing") of the probe components and contains artifacts. To 
reduce the probe-ringing time a low Q-factor NQR probe 
has been used with Q« 10-20. So, the signal can only be 
detected after a minimum time tjj ('dead time') which 
usually increases markedly with decreasing frequency. In 
our experiments the minimum practical r values varied 
from 70 [is at 1.5 MHz to 40 [is at 9 MHz. 



A. NQR Background 

If the crystal field symmetry is lower than cubic, 
the Hamiltonian of the quadrupolar interaction can be 
written^: 



n Q = ^(3Jl-/(/-i) + v (i x -4)), (i) 

where the quadrupole frequency 

= 3eQV ZZ 
VQ ~ 21(21 - l)h 

is defined by the nuclear quadrupole moment Q and the 
largest principal axis component Vzz of the EFG tensor 
and r[ = (Vxx — Vyy)/Vzz is the asymmetry parameter 
(here the principal axis of EFG tensor are selected such 
as \V Z z\ > \V Y y\ > \V X x\)- In NQR experiments the 
magnetic transitions between energy levels with Am = 
±1 are observed. 



B. Experimental techniques 

The NQR measurements were done using a home-built 
coherent pulsed NMR/NQR spectrometer. NQR spectra 
were taken "point by point" with a 7r/2 — r — ir radio 
frequency (RF) pulse sequence by varying the spectrom- 
eter frequency in equal frequency steps in the range 1.5- 
15 MHz at 4.2 K. The detailed NQR spectra were then 
constructed using a Fourier mapping algorith m 16 ' 17 . 



C. 2i Na NQR spectrum 

The nuclear spin energy levels splitted by the EFG are 
doubly degenerate (see Eq. [T|). So for 23 Na for which 
/ = 3/2 a single resonance line which corresponds to the 
±i <-> ±| transition can be observed at the frequency: 



v = v Q \jl + ^ 2 - (2) 

The NQR spectrum of 23 Na in the Na 2 / 3 Co02 com- 
pound measured at 4.2 K is shown in Fig. [3] It consists 
of three well resolved narrow lines which correspond to 
three Na sites with distinct local environment. Such a 
spectrum with resolved narrow lines undoubtedly con- 
firms that this phase is well ordered, as was pointed out 
in the earlier NMR study which has established that a 
characteristic signature of the x — 2/3 phased is to ex- 
hibit three distinct sodium sites. The shape of the 23 Na 
NQR lines has been fitted by a Lorentzian function and 
their positions and lincwidths (which were taken as a full 
width at half maximum - FWHM) are collected in Ta- 
ble El 

As for 7 = 3/2 the NQR frequency determined by 
Eq. ^ depends on both vq and ry, it is impossible to de- 
termine both parameters from the NQR spectrum. How- 
ever, these parameters were determined by NMR for each 
sodium site and are also listed in Table |TTJ Using these 
values and following Eq. © one can calculate the ex- 
pected frequencies f ca i c of 23 Na NQR lines (Table HI). As 
one can see the calculated and experimental values are 
in perfect agreement - therefore the notation proposed 
in Ref. Il4 for different sodium sites have been used in 
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FIG. 3: 23 Na NQR spectrum in the Na 2 / 3 Co02 compound 
measured at 4.2 K. The fully resolved quadrupole structure 
for the three Na sites is discussed in the text. 



Table [HI The tiny (<1.3%) difference in the experimen- 
tal and calculated frequencies for the Na2b site can be 
easily explained by a small deviation (not more than 5°) 
of the Z principal axis of the EFG tensor and the c crys- 
tallographic axis for this site. The relative intensities of 
sodium lines will be discussed later in section IV of this 
paper. 



D. 



3 Co NQR spectrum 



Former 59 Co NMR data*^ taken on this phase have 
already allowed us to evidence distinct Co NMR lines, 
but the NMR spectra were somewhat difficult to ana- 
lyze fully, as one needs to determine altogether the EFG 
parameters and NMR shifts of the various sites. The 
experiments were furthermore complicated by the need 
of a quasi perfect alignment of the powder sample with 
respect to the applied magnetic field. They however 
allowed us to determine some of the EFG parameters, 
which allowed us to know beforehand the range of ex- 
pected vq values. After a presentation of the spectra, we 
explain in some detail how the various transitions per- 
taining to the same sites can be identified. The identifi- 
cation of the Co NQR lines is explained in some details, 
and allow us then to determine the EFG values. 



1. Experimental results 

Typical NQR spectra of 59 Co in the Na 2 / 3 Co0 2 com- 
pound are shown in Fig. 2] We found that the number 
of lines and their intensities in the observed spectrum 
depend strongly on the delay r between pulses. This 
corresponds to the fact that the experimental spin-echo 
intensity depends on the rate of the spin-spin relaxation 
process. Consequently the fast relaxing nuclei are not 



TABLE II: Comparison of 23 Na NQR and NMR results. v caic 
is the calculated value for the NQR resonance line frequency 
using Eg. Q for the values of vq an n obtained by NMR (see 
text). 



Site 


NQR 
i/(MHz) Av (kHz) 


NMR 
vq (MHz) 


14 

n 


[/caic(MHz) 


Nal 

Na2a 

Na2b 


1.645(1) 8.2(1) 
1.928(2) 8.2(1) 
2.117(1) 9.0(1) 


1.645(5) 
1.74(1) 
1.86(1) 


0.01(1) 
0.84(2) 
0.89(2) 


1.645 
1.93 
2.09 



observable at long enough delay r between pulses, as 
will be detailed in section IV. In the spectrum observed 
with the shortest possible t$ = 45 [as there are 6 nar- 
row intense lines and several lines with weaker intensity 
(Fig. 2]). At the same time in the spectrum measured 
with tl — 100 lis the intensities of the high frequency 
lines greatly decrease, while the intensities of 4 low fre- 
quency narrow lines does not change a lot. So, it is ob- 
vious that in this phase with x — 2/3 two different types 
of Co sites - fast and slow relaxing - coexist, as was es- 
tablished already by NMRi£. 

As was shown earlier^i^ for the high sodium con- 
tent range (x > 0.5) of the sodium cobaltates phase 
diagram, a cobalt charge disproportionation into Co 3+ 
and Co + is a quite common picture. Cobalt ions in 
sodium cobaltates are in low spin configurations, so Co 3+ 
has an electronic spin S = and relax rather slowly 
in comparison with the cobalts with higher charge state 
on which holes delocalize (formally Co 4+ should have 
5=1/2). 



2. Identification of the NQR lines 

As the 59 Co nuclei have I — 7/2, for a single crys- 
tallographic site one should observe up to 3 lines in the 
NQR spectrum. The positions of these lines depends on 
both vq and r\. Figure [5^ shows the theoretical depen- 
dence of NQR frequencies of the allowed transitions on 
the asymmetry parameter of the EFG tensor for nuclear 
spin 7/2. These data were obtained by numerical diag- 
onalization of the main Hamiltonian recalled in Eq. |1| . 
One can see that the frequencies of the ±| ±| and 
±| <-» ±| transitions weakly depend on the asymmetry 
parameter 77, contrary to that of the ±4 <-> ±| transi- 
tion (Fig. [5^,). For 77 = 1 only two resonance lines in the 
spectrum should be observed. Measuring experimentally 
the resonance frequencies for different transitions should 
allow then to determine both vq and rj for / = 7/2. 

Since the number of lines in the experimental 59 Co 
NQR spectrum (Fig. [3]) is much larger than 3 (number 
of the observable NQR lines for / = 7/2), an analysis 
has to be done to determine the triplet of lines which are 
associated with a given cobalt site in the unit cell of this 
x = 2/3 compound. The basic data which can be used for 
such an analysis are the theoretical ratios of frequencies 



6 



co 



3 



-Q 

i— 

re 



co 
c 




o 

CD 

'q. 
C/5 



- I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [ 

Cola 



non-saturated 




saturated 
Co2b 



Co1b 
Co2a A 



3.5 4 3.6 



4.2 




JU i 



I T g 



Cola colb 
_J 



J I 

J L. 



Na 2/3 Co0 2 T=4.2K 



Co2a 



Co2b 

x =45p.s 



J 



x L =100nS 



5 6 
Frequency (MHz) 



FIG. 4: Main panel: 59 Co NQR spectra taken at 4.2 K for short (rs = 45 /j,s) and long (tl = 100 fj,s) delay between pulses. 

Insert: the two weak intensity lines of the ±3/2< > ±5/2 Co2a and Co2b sites are revealed by saturation of the high intensity 

slow relaxing lines. 



of different transitions for a single cobalt site which are 
restricted as established from Fig. [5p. 

Let us consider as an example the two intense lines at 
«7.52 MHz and sa6.47 MHz (Fig. 0} which corresponds 
to an an experimental frequency ratio sal. 16. Although 
we know from NMR experiments^ that those EFG values 
correspond to distinct Co sites, we want to show here that 
a simple analysis of the NQR spectra does confirm that 
point. In theory such a ratio could be found for a single 
site with EFG's with 77 = 0.47 or 77 = 0.85 (see Fig. Hp). 
As r) — 0.47 corresponds to the ratio of frequencies of 
±| <-> ±| and ±~ *-* ±| transitions, then one should 
find in the spectrum a line corresponding to the ±| <-> 
±| transition at a frequency w 7.52*1.578 w 11.87 MHz. 
We indeed didn't find such a line in our experiment. Also 
we know from NMR that such large EFG values do not 
exist^. Similarly for 77 = 0.85 one should find a line 
for the ±| <-> ±| transition at 7.52/1.51 w 4.98 MHz 
which does not exist either. Therefore the sa7.52 MHz 
and sa6.47 MHz lines cannot be attributed to the same 
cobalt site and correspond to two different fast relaxing 
cobalt sites Co2b and Co2a. 

Then one has to consider the next pair of lines, e.g. 
at «7.52 MHz and «4.83 MHz. The frequency ratio is 
equal sal. 56 and 3 possible values of 77 should be con- 
sidered. However only one value of 77 sa 0.36 gives the 
right 3rd line position in the spectrum, a weak and broad 
line at sa3.52 MHz which completes the spectrum of the 
Co2a site. We found this line by observing the slight 



asymmetry of the shape of the intense and narrow slow 
relaxing line at 3.581 MHz decreases with increasing r. 
To better resolve this line we used the fact that slow and 
fast-relaxing cobalts can be separated also by the large 
difference in their spin-lattice relaxation rates, as seen 
by NMR and as detailed later in section IV. To saturate 
the signal an additional 7r/2 pulse with 100/zs delay was 
used before the ir/2 — t — tt pulse sequence. In the spec- 
trum obtained after such a pulse sequence (" saturated" in 
the insert of Fig. 2]), the intensities of the slow-relaxing 
cobalt lines decrease and the lines of the fast-relaxing 
cobalts can be better resolved. This method allowed us 
to detect one more fast relaxing line at 4.145 MHz which 
completes the three line spectrum for the 2nd fast relax- 
ing Co2a site, with the two other lines at 6.473 MHz and 
3.03 MHz. 



Such an analysis allowed us also to determine that the 
two slow-relaxing lines at 3.581 and 4.178 MHz are the 
±| <-> ±| transitions of two different sites of cobalt, 
Cola and Colb, as could already be anticipated from the 
Co NMR data»i£ The corresponding ±| ±| transi- 
tions are at 2.387 and 2.785 MHz. We did not attempt to 
observe the ±i <-> ±| transitions for these slow-relaxing 
cobalts which should appear at low frequencies sa 1.19 
and sa 1.39 MHz, outside of the frequency range of our 
spectrometer. 
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TABLE III: Measured frequencies for the NQR transitions of 
the four Co sites in Na2/3CoC>2 and deduced values of vq and 
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rj - asymmetry parameter 



FIG. 5: (a) Theoretical dependence of the allowed NQR tran- 
sition frequencies on the asymmetry parameter rj of the EFG 
tensor for a nuclear spin / = 7/2. (b) Frequency ratios of 
the allowed transitions plotted versus rj. The particular case 
of the ratio 1.16 used in the text to illustrate the analysis is 
shown by dotted lines. 



' Co EFG parameters 



All the observed 59 Co NQR lines were fitted by 
Lorentzian functions and their positions and linewidths 
are collected in Table IIIII and Table IVI1 respectively. 
The values of vq and rj found in this work are in per- 
fect agreement with those quoted in the Ref. O. The 
fact that we observe a limited number of narrow and 
well-resolved 59 Co NQR lines (Fig. [4j) confirms the exis- 
tence of long-range order in this x — 2/3 phase of sodium 
cobaltates. Moreover NQR is clearly able to distinguish 
the charge segregation between the slow relaxing non- 
magnetic Co 3+ (sites Cola and Colb), and the fast re- 
laxing cobalt sites Co2a and Co2b, on which holes are 
delocalized. Despite all our efforts we have not found in 
the 59 Co NQR spectrum any traces of the Co3 site an- 
ticipated from NMR (see Ref. [l8l ). As explained in the 
Ref. [l] we did conclude that this apparent "site" was an 
experimental artifact due to imperfect alignment of the 
NMR sample. This point will be detailed in an extensive 
experimental report of the Co NMR spectral. The in- 
tensities of the 59 Co NQR lines will be discussed later in 
section IV after the analysis of the spin-spin relaxation 
of these cobalt NQR signals. 



Site 


Frequency ( MHz ) 
±i~±f ±§~±| ±f ~±| 


^q(MHz) 


n 


Cola 
Colb 
Co2a 
Co2b 


2.387(1) 3.581(2) 
2.785(1) 4.178(1) 
3.03(1) 4.145(3) 6.473(1) 
3.52(1) 4.826(2) 7.524(1) 


1.193(1) 
1.392(1) 
2.187(1) 
2.541(1) 


<0.017 
<0.016 
0.362(5) 
0.358(4) 



4- NQR linewidths 

The broadening of the Col and Co2 groups of lines 
were found to display different behaviours likely linked 
with their different magnetic properties. Figure O shows 
the linewidths variation with resonance frequency of the 
various transitions for each of the cobalt NQR lines. 
Two different behaviors appear clearly: the linewidths in- 
crease with increasing frequency for Col sites while the 
opposite holds for the Co2 sites. The Cola and Colb 
sites are non-magnetic with localized holes and one ex- 
pects essentially a broadening due to a distribution of 
EFG values. Indeed for a spin 7/2 system with zero 
asymmetry parameter (Col case), a spread eAq of EFG 
values induces broadenings of the different NQR transi- 
tions such that Avi/vi cx Aq/q, so that the linewidth 
increases linearly with the frequency of the transition as 
found in the data of Figure [5] One would expect then a 
linewidth of « 6 kHz for the lowest (±i <-> ±|) transi- 
tions of the Cola,b sites. Extrapolating the linear fits for 
Cola,b sites to zero frequency gives a very small residual 
value « 1 kHz, which could be attributed to dipole dipole 
coupling as will be detailed in section IV. 

As Co2a and Co2b are magnetic sites, a broadening 
due to magnetic interactions might take place and appar- 
ently increases with decreasing index of the quadrupole 
transition. This will be discussed in section IV after con- 
sidering the spin spin T2 processes. 



E. Rough estimates of the EFG values 

We have calculated the EFG tensors and the corre- 
sponding parameters of the quadrupole Hamiltonian (vq 
and rj) for all sites of the 23 Na and 59 Co in the simplest 
point charge model assuming charges 1 + on Na, 3 + on 
Col, 3.44+ on Co2 and 2" on O ( see Table [TV}. Of 
course one does not expect such point charge calcula- 
tions to fit perfectly the data as similar calculations do 
not even explain the simple case of NaiCo02 for which 
all Co are Co 3+ (Ref. lfj| and references therein). Such an 
approach does not take into account correctly the cobalt- 
oxygen covalency. We do however clearly see in the re- 
sults displayed in the Table QV] that for the Nal and Col 
sites the asymmetry parameter is of course that expected 
from the site symmetry, that is 77 = 0. We also find that 
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FIG. 6: Lines widths at half height for the NQR transitions 
of the different sites. The main panel is for Cola,b and inset 
for Co2a,b sites. 



the asymmetry of the structure yields rather large asym- 
metries for the EFG values on the Na2 and Co2 sites 
as found experimentally. However the magnitudes of the 
EFG computed do not agree quantitatively with the ex- 
perimental values. We did however find that the dis- 
placement of the Na atoms found in the X-ray structure 
analysis yields significant changes in the calculated EFG 
values. This reveals that more sophisticated calculations 
taking into account the atomic displacements and the 
band structure at least in the local density approxima- 
tion (LDA) are required to explain the EFG parameters. 



IV. SPIN-LATTICE RELAXATION AND 
SPIN-SPIN RELAXATION 



When thermal equilibrium of the nuclear spins is dis- 
turbed by RF pulses, their equilibrium magnetization is 
recovered by the nuclear nuclear spin-lattice relaxation 
(NSLR) process which corresponds to the relaxation of 
the longitudinal component of the magnetization. The 
decay of the transverse component of the magnetization, 
the nuclear spin-spin relaxation (NSSR) process is con- 
nected with the loss of coherence inside the spin system 
due to spin-spin interactions. Generally both relaxation 
processes could have different origins which reflect mag- 
netic and electronic properties of the materials, the inner 
arrangement, the interactions in the spin system, the dif- 
ferent motion and diffusion processes. In this section we 
report and discuss first in (A) the 59 Co spin-lattice re- 
laxation for all four cobalt sites. Next we consider in (B) 
the spin-spin relaxation starting from 23 Na as a simple 
case of NSSR. Then we analyze the 59 Co spin-spin re- 
laxation. We demonstrate that both NSLR and NSSR 
allow to differentiate the two types of Co sites - non- 
magnetic Co 3+ ions on 25% of the cobalt sites and a 
metallic Kagome network of Co sites on which the doped 
holes are delocalizedi 



TABLE IV: Results of the point charge model calculations 
of the EFG parameters for established structure of the 
Na 2 /3Co02- All quadrupolar frequencies vq are in MHz. 8 
(in degree) is the angle between the Z principal axis of the 
EFG tensor and the c crystallographic axis. Na: 23 Q=0.1; 
23 7=-4.1; Co: 59 Q=0.42; 59 j=-7; 



Site 


Experiment 


Non-shifted 


shifted 


Nal 


1/0=1-645(5) 
r) = 
= 


vq = 0.61 
77 = 
9 = 


Vq = 0.28 

7j = 

9 = 


Na2a 


i/ Q =1.74(l) 
77 = 0.84(2) 
9 = 


Vq = 0.93 
r) = 0.88 
9 = 90 


Vq = 0.94 

77 = 0.19 
9 = 90 


Na2b 


I/g=1.86(l) 
77 = 0.89(2) 
9 < 5 


vq = 0.53 
r) = 0.31 
9 = 90 


v Q = 1.15 
77 = 0.36 
9 = 90 


Cola 


v Q =l. 193(1) 
77 < 0.017 


vq = 1.06 
77 = 
9 = 


v Q = 1.71 
ij = 
9 = 


Colb 


i/q=1.392(1) 
77 < 0.016 


VQ = 1.18 
77 = 
9 = 


v Q = 0.94 
77 = 
9 = 


Co2a 


1-0=2.187(1) 
77 = 0.362(5) 


vq = 1.42 
f\ = 0.26 
9 = 18 


vq = 1.55 
r] = 0.56 
9 = 14 


Co2b 


i/q=2.541(1) 
77 = 0.358(4) 


vq = 1.55 
77 = 0.56 
9 = 11 


vq = 1.37 
77 = 0.40 
9 = 12 



A. Co spin - lattice relaxation 

To study the NSLR process we have used the mag- 
netization inversion recovery method which uses three 
pulses: 7r — t — 7r/2 — t — 7r, where t and r are the time 
intervals between pulses. In this sequence the first pulse 
rotates the magnetization by 180 degree, the second and 
third pulses gives a spin-echo, and the dependence of the 
spin-echo intensity on delay time t allows to monitor the 
recovery of the nuclear magnetization associated with a 
given NQR transition: 



M(t) = M (l-B-R(t)). 



(3) 



Here Mq is the thermal equilibrium value of magnetiza- 
tion and the parameter B characterizes the actual magne- 
tization after the first pulse at t — (the imperfection of 
the experimental conditions gave typical values B ~ 1.8 
rather than B = 2 expected for a perfect it pulse). The 
shape of the relaxation function R(t) depends on the nu- 
clear transition sampled. For a two-level nuclear system 
(like the 7 = 3/2 NQR case) this process is exponential 
R(t) — exp(—t/T\), with a characteristic time constant 
Ti, the spin-lattice relaxation time^i. But generally for 
I > 1/2 the nuclear energy levels are differentiated by 
the quadrupole interaction with the crystalline electric 
field (Eq. [T|) . As a consequence the difference in popula- 
tion between adjacent levels which are probed by the RF 
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TABLE V: Theoretical spin-lattice relaxation function Q pa- 
rameters for different 59 Co sites for the case of magnetic re- 
laxation by weak isotropic fluctuating magnetic fields (based 
on Ref.[M). 



1 F - ' r- 



Site 


V 


Transition 




Ai 


2 


A 2 


«3 


A 3 


Cof 









0.14 


21 


0.65 


10 


0.21 


3 






0.74 


21 


0.16 


10 


0.1 


3 


Co2 


0.36 


±1- 


0.27 


17.7 


0.53 


8.96 


0.2 


3 








0.88 


17.7 


0.038 


8.96 


0.082 


3 



pulses depends on the populations of the levels which are 
not hit by the RF pulses. Therefore the magnetization 
recovery becomes multi-exponential 
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R(t) = E 



a.texp 



(4) FIG. 7: Spin-lattice relaxation curves for Co2b site measured 



but is still characterized by a single T± value. 

The spin-lattice relaxation could be driven either 
by magnetic or quadrupolar fluctuations. However in 
the sodium cobaltates, which exhibits unusual magnetic 
properties, the magnetic relaxation mechanism domi- 
nates at least at low T. Following Ref. Hi] we obtained 
the parameters for the theoretical relaxation functions 
R(t) for the case of magnetic relaxation by weak isotropic 
fluctuating magnetic fields for different transitions of dif- 
ferent Co sites - see Table El 

Figure [7] shows experimental spin-lattice relaxation 
curves for the ±| *-> ±| and ±| «-> ±| transitions for 
the Co2b site measured at 7.52 and 4.83 MHz, respec- 
tively. The fits with Eq. (4) with the ai and Ai taken 
from Table El are shown by solid lines. These fits give 
T\ values of 0.73(4) ms and 0.8(1) ms respectively for 
the ±| ±| and ±| ±| (Table ED). This rather 
good agreement confirms that the relaxation is caused 
by magnetic fluctuations and that the spectral density of 
fluctuating fields is nearly independent on the frequency. 
For the Cola and Colb sites it was found that T\ values 
±| is shorter than that for ±| <-> ±| 



for the ±| 



(Table fVTj) . This effect does result in that case from the 
proximity of the fast relaxing Co2 ±| <-> ±| transitions 
(see Fig. [4| which induce by cross-relaxation processes 



a reduction of the the recovery time on the ±| 



transitions of Cola and Colb. For these Col sites the Xi 
value retained then hereafter is the longer one deduced 
from the ±| •*->• ±| transition data. 

In Fig. [5]we show the temperature dependence of nu- 
clear spin-lattice relaxation rates 1 jT\ for all four cobalt 
sites. As one can see for both Co2a and Co2b NSLR is 
about two orders of magnitude larger than for Cola and 
Colb. As was already stated in Ref. [l|, this proves the 
non-magnetic nature of Col sites. Also Fig.[5]shows that 
1/Ti increases with T at low temperatures. Above 20 K 
the NSLR becomes so fast that we loose then the Co2 
NQR signals. Measurements at higher T cannot be car- 
ried out then by NQR, and T\ evolution will be rather 
studied by NMR as will be discussed in a forthcoming 



at 7.52 MHz and 4.83 MHz (±| <-> ±| and ±| <-> ±| transi- 
tions correspondingly) at 4.2 K. Solid lines are the fits of ex- 
perimental points by the magnetization relaxation functions 
@ with parameters from the Table IVl 



paper—. Similarly the spin-lattice relaxation measure- 
ments for 23 Na were technically challenging and have 
therefore not been performed here as they have been 
studied in great detail by NMR in Ref. |JJJ for various 
Na x Co02 phases including the x = 2/3 phase. 



B. Spin-spin relaxation and signal intensities 

NSSR is a complex phenomenon, but at its most fun- 
damental level the random fluctuations of the local mag- 
netic field leads to a loss of the initial phase coherence 



0.08 
0.06 
0.04 
0.02 



±5/2-±7/2 
o Co2a 
□ Co2b 



±3/2-±5/2 
a Cola 

v Colb 







8 10 12 14 

Temperature (K) 



20 22 



FIG. 8: Temperature dependence of nuclear spin-lattice re- 
laxation rates 1/Ti for all four cobalt sites in the Na 2 /3Co02 
compound. For Cola and Colb measurements were done on 
the ±| <-» ±| transitions and for Co2a and Co2b measure- 
ments were done on the ±| «-> ±| transitions. Linear fits 
shown by dotted lines are guides to the eyes. 
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TABLE VI: 59 Co relaxation and NQR linewidth parameters (see text for details). 



Site Cola Colb Cola Colb Co2b Co2a Co2b 

Transition ±f <-> ±f ±jj-> ±f ±| <-> ±| ±j g ±j ±| g ±j ±§ <-> ±| ±j <-> ±| 

j/ (MHz) 2.387(1) 2.785(1) 3.581(2) 4.178(1) 4.826(2) 6.473(1) 7.524(1) 

Ti (ms) 40.5(5) 55.5(5) 26(1) 46(2) 0.8(1) 0.59(3) 0.73(4) 

T 2 (us) 316(4) 326(7) 320(2) 335(6) 30(2) 71(2) 49(2) 

n 2 2 2 2 1 1 1 

T{ (ms) 2.3 3.2 2.6 4.6 0.050 0.058 0.072 

a 0.61(5) 1.20(8) 0.67(5) 

Av ca i c (kHz) 1.0(4) 0.9(2) 1.0(2) 0.9(2) 10.6(3) 4.5(3) 6.5(2) 

Ai/(kHz) 12.8(9) 12.1(8) 18.6(9) 18.3(8) 42(2) 24(1) 24(1) 



of the nuclear spins and therefore to the decrease of the 
transverse nuclear spin magnetization. In the Redfield 
theory^ the transverse relaxation arises from two mecha- 
nisms: (1) nuclear spin-spin coupling (via magnetic dipo- 
lar or indirect interactions) and (2) the spin-lattice relax- 
ation. 

The NSSR is studied by monitoring the spin-echo in- 
tensity as a function of the time delay r between ir/2 and 
7r pulses. In general the decay of transverse magnetiza- 
tion M(t) as a function of time t = 2r can be fitted by 
the equation: 



.V/(7) = .U(())r, 7 ,(-(± 



(5) 



The NSSR process is characterized by the relaxation time 
T2 and the exponent n, which characterizes the decay 
shape which usually varies between Lorentzian (n = 1) 
and Gaussian (n — 2). 



1. 



Na spin-spin relaxation 



The spin-spin relaxation decays for sodium sites are 
shown in Fig. [9l The shortest r value used was only 
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FIG. 9: Spin-spin relaxation curves of 
the Eq. ©. 



Na sites and fits by 



~ 65 fj,s due to the long 'dead time' Td of the spectrome- 
ter (see section III) . Then the beginning of the relaxation 
decays were lost and only tails were measured. There- 
fore it was impossible to determine reliably the transverse 
magnetization decay shape. Nevertheless the experimen- 
tal data were fitted by a Lorentzian (n = 1) function (see 
Eq. iJSJ)) in order to obtain estimates of the spin-spin re- 
laxation times and of the NQR lines intensities. The T\ 
values of 23 Na known from the data of Ref. [lfj for the 
same phase of sodium cobaltates were long enough to as- 
certain that they do not contribute significantly to the 
transverse magnetization decay. So the values obtained 
by fitting the relaxation decays for the sodium sites given 
in Table [VTTI can be associated with spin-spin processes. 

The values of magnetization at zero time M(0) after 
v 2 frequency correction allowed us to estimate the rela- 
tive intensities of sodium lines (sites) which are also re- 
ported in Table EE These data are in good agreement 
with the more accurate ones obtained by NMR. Here the 
accuracy^ is indeed limited primarily by the extrapola- 
tion required from Try to r = 0, and also by the variation 
of spectrometer sensitivity with frequency. 

Due to the NSSR process the linewidth of an NQR sig- 
nal at half-height should be equal A^ ca ; c = l/(7r • T2^- 
Such Av ca i c for different Na sites are also collected in 
Table IVIII And it is easy to see that they are approxi- 
mately three times smaller than the experimentally mea- 
sured linewidths for all sodium sites (Table IVIip . This 
clearly shows that the experimentally observable Na lines 
are inhomogeneously broadened by a weak distribution 
of EFG such that Avq /vq rs 4 • 10~ 3 . Such a small value 
found indeed imply a very well ordered structure in the 
sodium layers of this Na 2 /3Co02 compound. 



TABLE VII: 23 Na relaxation and NQR linewidth parameters. 



Site 


/(%) 


T 2 (jjs) 


k-Vcalc (kHz) 


Au (kHz) 


Nal 


23(5) 


125(2) 


2.5(1) 


8.2(1) 


Na2a 


46(8) 


82.3(7) 


3.8(2) 


8.21(5) 


Na2b 


31(7) 


97(1) 


3.3(2) 


9.02(6) 
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FIG. 10: Spin-spin relaxation curves for Cola, b sites (a) on 
±3/2 <-> ±5/2 transition, and for Co2a, b (b) on ±5/2 <-> 
±7/2 transition. 



2. 



' Co spin-spin relaxation 



The transverse relaxation decays for 59 Co sites are rep- 
resented on Figure[l0l where panel (a) shows data for the 
±| •*->• ±| transitions of the Cola,b sites and panel (b) 
that for the ±7/2 <-* ±5/2 transitions of the Co2a,b sites. 
The results of the fit of these curves by Eq.[5]are summa- 
rized in the Table I VII For the Col sites the transverse 
magnetization decay has a Gaussian form [n — 2), for 
Co2 it is a Lorentzian with n = 1. To understand this 
difference one needs to analyze the possible influence of 
spin-lattice relaxation on the transverse relaxation. 

For the cobalts the spin-lattice relaxation is fast 
enough so its contribution to the spin-spin relaxation be- 
comes important. As the spin-lattice relaxation of 59 Co 
is multi-exponential (Eq. 0]) then a good approximation 
might be done by considering that Ti in the Eq. [5] should 
be related to the slope of the NSLR curves T[ at short 
times. For example, the Eq. [4] with the parameters from 
the Table [V] can be written, at short enough times for 
the Co2b's ±| <-> ±|transition as 



B{t) 



10.12* 



Ti 



= 1 - 



t 



In the Table I VII one can see the estimated values of T[ 
for different Co sites and transitions. For Cola and Colb 
these T[ times are much longer than the spin-spin relax- 
ation Ti on Fig.fTUk. Therefore the NSLR contribution to 



the NSSR is negligible and we see in experiments the real 
spin-spin relaxation of Col sites. But for Co2, the calcu- 
lated T[ values become comparable with the obtained Ti 
times. To emphasize this we reported in the Table lVTl thc 
ratio a — TijT' x . As one can see a is almost the same for 
the two measured transitions ±| <-»• ±1 and ±| <-> ±| 
for the Co2b site. So we conclude that for Co2 sites the 
measured Ti values are determined by spin-lattice relax- 
ation process. This explains as well the difference in the 
transverse magnetization decay shape for Col and Co2 
sites. 

This influence of spin-lattice relaxation on the spin- 
spin relaxation explains why the intensity of the low fre- 
quency transitions of Co2a and Co2b sites are found 
so weak experimentally (see Fig. [4]). As one can see 
from the Table IVIl for the Co2b site, the measured 
T2 =30(2) us for the ±| «-> ±| transition is smaller 
than that (49(2) /us) for the ±| <-> ±| transition. Also 
as was already stated in section III, the minimum us- 
able time td between RF pulses increases with decreas- 
ing frequency. So both factors reduce the intensities of 
the observed low frequency Co2 NQR lines. 

As has been noticed previously in section III, the sig- 
nal intensities observed at different frequencies are af- 
fected by the variation of spectrometer sensitivity, so 
careful determination of the relative site occupancies of 
Col and Co2 could not be done directly from our ex- 
periments. But we could estimate the fractional occu- 
pancies inside each group of cobalt (Col and Co2) sep- 
arately, as their resonance frequencies are close. After 
correcting for the spin-spin relaxation T2 and for the 
v 2 frequency dependence of the signal intensity we ob- 
tained: Cola/Colb=l. 95(0.1) and Co2b/Co2a=1.9(0.2), 
in agreement with similar determinations from 59 Co 
NMR, which allowed as well better determinations of the 
Col/Co2 ration. 

Spectral lines of cobalt are inhomogeneously broad- 
ened (see Table IVTl) and this is clear if one compares the 
experimental linewidth with the calculated one. Notice 
that the ratio Ai//Az/ ca/c is w 4.0(0.4) for both ±| «-► ±| 
and ±| «-> ±|) transitions of Co2b and w 5.0(0.3) for 
Co2a, which confirms that these broadenings have the 
same origin. 



V. CONCLUSION 

In summary, we have reported detailed powder X ray 
and NQR data for the Na2/3Co02 compound support- 
ing the recent paper— in which we proposed a detailed 
structural model for this material from NMR/NQR re- 
sults. The NQR spectrum of 23 Na and 59 Co consists of 
several narrow lines, which could be associated with 3 
sodium and 4 cobalt sites, in a 2D unit cell comprising 
8 Na sites over 12 Co. The finite number of unequiva- 
lent site positions confirms that the simple 3D ordering of 
the Na layers leads to a differentiation of cobalt sites. We 
have performed as well a careful Rietveld analysis of the 
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X ray data which confirms the atomic 3D unit cell pro- 
posed from NMR/NQR. The best Rietveld fits allowed us 
to determine the slight atomic displacements of the Na 
atoms in the model unit cell. Also the quality of the fit 
allows us to ascertain the x=2/3 value of the Na content. 

Others have also detected this phase and identified it 
with the present one from its large spin susceptibility at 
low temperatures 11 . They have performed single crys- 
tal diffraction experiments in which they also find a 12 
Co 2D unit cell. Using chemical analyzes of their single 
crystals they proposed that x — 0.71, and have elabo- 
rated then a structural model in which the concentration 
of Na alternates between planes with x = 2/3 with 8 
Na per 2D cell and x = 3/4, with 9 Na per 2D cell. 
In these Na layers the Na vacancies would be alterna- 
tively organized in divacancies and trivacancies. One can 
immediately see that such a structure would never ex- 
plain the present NMR/NQR data as it corresponds to 
a larger number of distinct Na sites (3 in each Na plane) 
than found experimentally. Furthermore the only sites 
exhibiting the threefold symmetry and r\ = would then 
be the Nal sites of the x = 2/3 planes. The fraction of 
such sites would be x = 2/17, nearly twice smaller than 
the 25% = 2/8 found experimentally. 

The spin-spin and spin-lattice relaxation of the NQR 
lines were studied for all Co sites, and the data analysis 



has allowed us to demonstrate that both T\ and T% data 
do allow to establish the difference of magnetic properties 
of the two types of Co sites. As indicated in Ref. [0, they 
constitute two sublattices: a Kagome structure for the 
magnetic sites and the complementary triangular lattice 
of non-magnetic Co sites involving 25% of the Co. 

Beyond the simple results obtained for this Na 2 /3Co02 
compound, the approach developed here will be cer- 
tainly useful to perform structural determinations for 
pure phases with different Na contents. We have in- 
deed isolated some of them, which display quite distinct 
ground state physical properties^. 



VI. ACKNOWLEDGMENTS 

We thank N. Blanchard for her help in synthesiz- 
ing materials and Y. S. Meng for numeourous exchange 
about their computations. We acknowledge financial sup- 
port by the ANR (NT05-441913) in France and by the 
grant RNP-6183 in Russia. Expenses in Orsay for I.R.M. 
and A.V.D. have been supported by the Triangle de la 
Physique. T.A.P. has obtained a fellowship from the E.U. 
Marie Curie program "Emergcntcondmatphys" for part of 
her PhD work performed in Orsay. 



* Electronic address: Irek.Mukhamcdshin@ksu.ru 

1 H. Alloul, I. R. Mukhamedshin, T. A. Platova, and 
A. V. Dooglav, EPL 85, 47006 (2009). 

2 I. Terasaki, Y.Sasago and K.Uchinokura, Phys. Rev. B 56, 
R12685-R12685 (1997) 

3 K. Takada, H.Sakurai, E.Takayama-Muromachi, F.Izumi, 
R.A.Dilanian and S. Sasaki, Nature (London) 422, 53 
(2003) 

4 G. Lang, J. Bobroff, H. Alloul, G. Collin, and N. Blan- 
chard, Phys. Rev. B 78, 155116 (2008) 

5 G. Lang, J. Bobroff, H. Alloul, P. Mendels, N. Blanchard, 
and G. Collin, Phys. Rev. B 72, 094404 (2005) 

6 C. de Vaulx, M.-H. Julien, C. Berthier, M. Horvatic, P. 
Bordet, V. Simonet, D. P. Chen, and C. T. Lin, Phys. 
Rev. Lett. 95, 186405 (2005) 

7 H. W. Zandbergen, ML. Foo, Q. Xu, V. Kumar, and 
R. J. Cava, Phys. Rev. B 70, 024101 (2004). 

8 P.H.Zhang, R.B. Capaz, M.L.Cohen, and S.G. Louie, Phys. 
Rev B 71, 153102 024101 (2005) 

9 J. Bobroff, G. Lang, H. Alloul, N. Blanchard, and G. 
Collin, Phys. Rev. Lett. 96, 107201 (2006) 

10 H. Alloul, I. R. Mukhamedshin, G. Collin, and N. Blan- 
chard, EPL 82, 17002 (2008). 

11 F. C. Chou, M.-W. Chu, G. J. Shu, F.-T. Huang, Woei 
Wu Pai, H. S. Sheu, and P. A. Lee, Phys. Rev. Lett. 101, 
127404 (2008). 

12 J. D. Jorgensen, M. Avdeev, D. G. Hinks, J. C. Burley, 
and S. Short, Phys. Rev. B 68, 214517 (2003). 



13 Q. Huang, M. L. Foo, R. A. Pascal, Jr., J. W. Lynn, 
B. H. Toby, Tao He, H. W. Zandbergen, and R. J. Cava, 
Phys. Rev. B 70, 184110 (2004). 

14 I. R. Mukhamedshin, H. Alloul, G. Collin, and N. Blan- 
chard, Phys. Rev. Lett. 93, 167601 (2004). 

15 A. Abragam, The Principles of nuclear magnetism, Ox- 
ford: Clarendon Press, London (1961). 

16 W. G. Clark, M. E. Hanson, F. Lefloch, and P. Segransan, 
Rev. Sci. Instrum. 66, 2453 (1995). 

17 A. P. Bussandri, and M. J. Zuriaga, J. Magn. Reson. 131, 
224 (1998). 

18 I. R. Mukhamedshin, H. Alloul, G. Collin, and N. Blan- 
chard, Phys. Rev. Lett. 94, 247602 (2005). 

19 I. R. Mukhamedshin, H. Alloul, in preparation. 

20 C. P. Slichter, Principles of Magnetic Resonance, 2nd ed., 
Springer, New York, (1980). 

21 J. Chepin and J. H. Ross, J. Phys.: Condens. Matter 3, 
8103 (1991). 

22 Y. Hinuma, Y. S. Meng, and G. Ceder, Phys. Rev. B 77, 
224111 (2008). 

23 We did not anticipate initially such a large 2D unit cell 
initially proposed by Hinuma et alw^, and searched a 
twice smaller commensurate lattice structure with a ort = 
ahe-x^fi, b or t = Sdhex, Cort = 3ch e ^- Although in this 
model one could index the most significant satellite lines 
correctly, the fit of intensities was not as satisfactory as 
with the complete substructure used. 



